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1. Introduction 

The study of nuclear reactions has prompted the delineation of fundamental princi- 
ples whose importance extends beyond nuclear physics, from where they were originally 
inferred, but reaches vast sectors of physics in general and perhaps of other fields as well. 

Nuclear reactions are of central interest because they deal with strong interactions, 
for which perturbative approaches are of little value. Indeed, these strong forces, acting 
at short distances between the nucleons, entail large matrix elements to excited states of 
the nucleus lying at high energy: this occurrence in turn leads to a poor convergence of 
the diagrammatic expansion. In reaction theory the above problem finds a solution in the 
recognition that once an average value (to be appropriately defined) is subtracted from 
these matrix elements, then the resulting quantities are random. Actually, we believe 
these results to hold not only for the continuum, but for the nuclear bound states as well 
and also for other systems beside the nuclear one. In other words we conjecture that 
many strongly interacting systems in states of high excitations are chaotic. 

In this conceptual framework an approach is here presented aiming at the derivation 
of the nuclear mean field (MF), binding the nucleons together in the nuclear ground 
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state, from statistical considerations. In general, the mean field happens to be a good 
description for a variety of many-body systems, whose constituents can interact either 
with strong or weak forces, of long or short range. It is thus natural to search for 
the features, common to the various systems, ensuring the validity of the mean field 
approximation. In line with the above speculations we hold the view that the statistical 
aspect is the common feature linking all these different systems together. 
Accordingly, we base our treatment of the MF on two propositions: 

a) the MF can be obtained through an energy average that smooths out short time 
events as, e. g., those corresponding to the violent short-range nucleon-nucleon 
encounters in nuclei; 

b) the matrix elements of the residual interactions, namely the force left out after the 
average interaction has been removed, are random with zero average value. 

As we shall see these two elements suffice to obtain explicit expressions for both the 
MF and the fluctuations (the error) away from the average. 

Clearly the error should be small for the MF to be meaningful. In this connection 
it turns out that, because of item b), the average error actually vanishes. However, the 
average of the square of the error is not vanishing and its magnitude we wish to assess. For 
this purpose, we expand the fluctuations of the MF in terms of excitations of increasing 
complexities, each class of excitations providing a contribution to the error. Notably, a 
general condition for the convergence of this expansion (which is a finite one in the case 
of atomic nuclei) can be derived; in fact, a parameter can be identified determining the 
relative importance of each contribution to the error. 

It is finally worth emphasizing that the theory here developed is, in the one hand, 
quite general in the sense that no assumptions on the nature of the system being studied 
are made, and, on the other, closely inspired from, in fact almost identical with, the 
statistical theory of nuclear reactions. Indeed it provides, as the latter does, an alternative 
method to account for the contributions to physical observables stemming from processes 
slowly and rapidly varying with the energy. It differs from the customary approaches 
which are based on the explicit consideration of the Pauli and dynamical correlations, 
the latter typically providing rapidly varying contributions. 

In order to get an appreciation on how the method is practically implemented an 
exploratory calculation of the ground state energy of an infinite, interacting Fermi gas, 
i. e. nuclear matter, will be discussed at the end of this lecture. For this system a 
reasonable estimate of the error is obtained. Furthermore, and significantly, the residual 
interaction appears to be drastically reduced with respect to the bare one, an occurrence 
directly related to the statistical nature of the present approach. 

The present lecture, which is based on Rcf. [1], is organized as follows. In sect. 2 the 
formalism is revisited, — in particular the partition of the Hilbert space into the P and Q 
sectors, — and the related equations are discussed, while in sect. 3 a precise definition of 
the energy average is provided. The latter is performed with a distribution, characterized 
by a single parameter e, amounting to shift the position of the poles corresponding to 
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the excited states of the system away from the ground state energy. In the same section 
we pave the way for setting up a suitable formalism to account for the corrections (the 
error) to the mean field energy. 

In sect. 4 the expression for the corrections to the mean field energy Eq is given in the 
framework of the expansion, above referred to, built out of successive 2p-2h excitations. 

In sect. 5 we deal with the issue of introducing an operator, projecting into a smooth 
P-space, simple enough to allow the present approach to be worked out for nuclear 
matter. We thus restrict the P-space, — also with the aim of rendering transparent the 
comparison between the mean field energy as given by the Hartree-Fock (HF) approach 
and by the present theory, — just to a simple state: The HF determinant. With this 
choice for P, we are able to establish an expression connecting the true ground state E, 
the mean field E and the Hartree-Fock Phf energies. Importantly, £0 is shown to be 
always lower than Phf and it turns out to be proportional to the variance, to be later 
defined, of the residual effective interaction V among nucleons, which, in our context, is 
precisely defined. 

In sect. 6 we address the issue of the occupation probability S 2 (spectroscopic factor) 
of the ground state of nuclear matter as defined by the structure of our P-space. We are 
able to deduce, within our scheme, an expression for S 2 in terms of E, Eq and Phf, and 
moreover embodying the energy derivative of V. 

Finally, in sect. 7 a scheme which permits numerical predictions of the present 
approach is developed. It refers to the infinite Fermi gas, namely nuclear matter. 

For sake both of simplicity and of illustration, we account for the first corrections 
only and, even so, a free parameter a should be introduced to obtain an estimate of their 
size. Two systems of two equations each can then be set up by coupling the equation for 
the mean field energy to the one for the fluctuations. By solving the systems and with a 
convenient choice of the two parameters entering into the present approach, namely e and 
a, results consistent with the empirical features of nuclear matter, as derived through 
an extrapolation from finite nuclei, are obtained. At the same time an estimate for the 
fluctuations of the mean field energy is provided. These goals are achieved with a residual 
effective interaction which, as previously emphasized, is lowered by about three orders of 
magnitude with respect to the bare one and yields a spectroscopic factor not in conflict 
with other independent estimates [3]. 

In the concluding section the approach is summarized and possible improvements 
discussed. 

2. Formalism 

In this section we shortly present the formalism originally introduced in Ref . [4] . Let 
H be the nuclear hamiltonian entering into the Schrodingcr equation 

(2.1) P* = P*. 

Let P be the hermitian operator projecting the nuclear ground state into the Hilbert 
subspace of functions associated with the low momentum transfer physics and Q the 
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operator complementing P. Clearly 

(2.2) P 2 = P, Q 2 = Q, PQ = 0, P + Q = l. 
It is then proved that the pair of equations 

(2.3) {E-Hpp)(P^) = Hp Q {Q^), 
and 

(2.4) (E - H QQ )(QV) = H QP (PV), 

are equivalent to (2.1). In the above formulas the following shorthand notations have 
been introduced, viz. 

(2.5) Hp P = PHP, H QQ = QHQ, H PQ = PHQ, H QP = QHP. 
Now a direct way to proceed in order to obtain PIS? is to get from (2.3) 

(2-6) (P*) = - 1 Hp Q (Q^) 

rL — tlpp 

and from (2.4) 

(2.7) (Qtt) = H QP( P *)- 
Inserting then (2.7) in (2.6) one obtains 

(2.8) {m = -J—H PQw J s -H QP m) 
or, equivalently, after multiplication by E — Hpp, 

(2-9) (e - Hpp - Hpq 1 H Q p) (Ptt) = 0, 

namely the equation obeyed by the component of \& in the P-space. A more general 
procedure is however required should a solution of the equation 

(2.10) (E-Hpp)^o = 

exists. 

In such a case, instead of eqs. (2.6) and (2.7), one has 
(2.H) (Ptt) = tt + - 1 H PQ (QV) 

Hi — / / pp 

and 

(2.12) (Q*) = - 1 H QP *o = —Hqp^o, 

E - Hqq - Wqq e Q 
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where 

(2-13) W QQ = Hqp 1 H PQ 

Hj — tl-PP 

and 

(2-14) cq = E — Hqq — Wqq . 

The above equations allow then one to recast (2.3) as follows 

(2.15) (E - Hpp){PV) = Hp Q —H Q p^o 

e Q 

and to express according to 

(2.16) * = (P*). 

1 + ~ Hp Q —H Q p 

E - Hpp eQ 

The combination of the two above equations leads in turn to 

(2.17) (E - Hp P ){P^) = Hp Q ^-H Q p (P*), 

l + E—lTp-p HpQ - Q HQP 

where no trace is left of ^o- Now using the operator identity 

(2 - 18) b TTcab = TTbca b > 

with the identifications 

(2-19) A = — , B = Hqp, C = 1 H PQ , 

eg E - Hpp 

one obtains 

(E - H P p)(PV) = Hpq-^— l - -H QP (PV) 

1 + Hqp— — H PQ — 

E - Hpp e Q 

(2.20) =H PQ - -H QP (PV). 

eQ l + W QQ ±- 

e Q 

The obvious identity 
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leads finally to 



(2.22) 



E-H 



pp 



H 



PQ- 



-H, 



V 



1 

eg 



QP 



— I +W QQ 



(P*) = 0, 



/ 



which is entirely equivalent to (2.9). 

It is important to note that the equation obeyed by P^f, being associated with the 
intricate many-body operator 



(2.23) 



H = Hpp + Hpq 



1 



-H QP , 



+ Wqq 



is not an eigenvalue equation in the usual sense. Indeed, the dependence upon E, the 
exact ground state energy of the system, is non-linear, since the latter appears also in 
the propagator in the Q-space in (2.22). 



3. Averaging upon the energy 

The partition of the Hilbert space into a P and a Q sector should be performed in 
conformity to the principle of including into the P-space the wave functions with the 
simplest structure, namely those with a smooth spatial dependence, and in the Q-space 
the wave functions of greater intricacy. The actual \& of the system should of course be 
viewed as a linear superposition of components with all the possible degrees of complexity. 

The problem then is: how to account for the average impact on the P-space of the 
Q-space wave functions, ignoring the detailed behaviour of the latter? To deal with this 
question an averaging procedure should first be prescribed. For this purpose we treat 
the energy E as a variable upon which both the wave function and the matrix elements 
depend. Clearly, the components of the wave function varying most rapidly with E lie 
in the Q-space: thus performing an average over the energy amounts to smoothing out 
the behaviour of Qty which can then be taken into account in the determination of the 
mean field. The latter rules the gentle physics taking place in the P-space. 

However, by replacing (Q'f) with (Q^f) we also change at the same time (PS&) into 
(P5») and E into, say, E (the angular brackets meaning, of course, energy averaging). 

The system equivalent to the Schroedinger equation will now accordingly read 

(3-1) (P^) = § + 1 H PQ (QM>), 

-C'O — tipp 

where 



(3.2) 



(Eq - H P p) * = 
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and 

(3.3) (Q*> = <— )H QP *o- 

e Q 

Inserting (3.3) into (3.1) yields 

(3.4) (Eo - Hpp) <P*) = H pq ( — )H Q p^o- 

Following then the same steps as before the equation 
/ \ 



(3.5) 



Eo - Hpp - Hpq — Hqp 

(— > + w QQ 

V e Q J 



(Ptf ) = (Eq - H) (P*) = 



is obtained: it clearly corresponds to equation (2.22) when a suitable energy average has 
been performed. 

We should now face the problem of specifying the energy average. For this purpose 
a smoothing function p(E, E ) with the property 

(3.6) J p(E,E )dE = 1 

should be introduced. Then averaging a function f(E) means to perform the integral 

(3.7) (/)= J p(E,E )f(E)dE 

which must be a real quantity since we are dealing with bound states. A smoothing 
function obeying the above conditions is 

^ ^ = ^-(le) ' 

In this case the integration in eq. (3.7) is performed along a path coinciding with the real 
axis ReE with a small semi-circle described positively about the singularity (Eq — e). 
If f(E) is bounded at infinity sufficiently strongly, the Cauchy's integral formula can be 
applied so that the condition of eq. (3.6) is satisfied and eq. (3.7) becomes 

(3.9) (f) = f(E -e). 

In the present approach e is an empirical parameter essentially measuring the range over 
which the average is taken. In other words, for a sufficiently large value the energy of 
the states within e can be neglected and it is in this sense that an average is performed. 

At this point the system corresponding to the original Schroedinger equation can be 
cast into the form (see also Kawai, Kerman and McVoy [5]) 

(3.10) (E-fL)(P^) = V PQ (Q^) 

(3.11) (E-H QQ )(Q*) = V QP (P^) 
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where 

(3.12) TL = Hpp + VpqVqp 
and 

(3-13) Vpq = H P Q^ 



Eq — e — E 



I Eo-e-E 
E Q - e - Hqq ' 



(3-14) v Q p = \ U E ° T h H Q p - 

y E - e - Hqq 

Note that eqs. (3.10) and (3.11), while having the same structure as the original eqs. (2.3) 
and (2.4), display a potential, coupling (P 1 ^) and (Q^), which is Vpq rather than Hpq. 
The strength of the coupling is thus considerably reduced by roughly (e/HQQ)^ which 
is much less than one. This in turn entails that H PP and TL are of the same order of 
magnitude. 

Now the spectral decomposition of the operator 1/(E — TL) can be performed in 
terms of eigenfunctions of the hermitian operator TL, namely 

(3.i5) m n = EA. 

It reads 

(Pi) = ^ WW = £ WW 
,3. 16 ) H4o) <5|^ + (_±_)V PQ < W . 

the prime on (1/(E — TL)) signifying that the lowest eigenfunction <J> is to be excluded. 
From eq. (3.16) it follows 

(3.17) <*o|P*> = <*f^*> . 
To obtain we return to eqs. (3.10) and (3.11) and get 

(3.18) (Qtf) = E _\ V Q p\Q ){$ \PV), 
where 

(3.19) h QQ = Hqq + Wqq 

and 

(3-20) Wqq EE V QP Vpq. 



A STATISTICAL APPROACH TO THE THEORY OF THE MEAN FIELD 



9 



Inserting eq. (3.18) into eq. (3.17) one finally obtains 

(3.21) E - E Q = (<5> Q \V PQ 1 Vqp\^o), 

E - riQQ 

where the right hand side represents the correction to the mean field energy Eq. 

Remarkably, the quantity (<^q\P^) does not appear in eq. (3.21), the reason being 
related to the arbitrariness of the amplitude of <&o. However, the bound state wave 
function <j> can be normalized: our choice is {^q\P^) = 1. 



4. Corrections to the energy of the mean field 

We now search for an expression yielding the corrections to the ground state energy 
£0 of the shell model in the statistical framework. This is accomplished along the same 
line adopted in the theory of nuclear reactions [2, 6]. Accordingly (3.21) is expressed 
through an expansion where each contribution is not identified by a coupling constant 
to some power, but rather by the degree of complexity of the states of the Q-space. 

One thus writes 



(4.1) AE = ^2 A^. 

where 



m—1 



(4.2) AE m = ($o\VoiG 1 V 1 2G2---V m - 1 , m G m V m , m - 1 G m - 1 ---V 10 \*o), 

(4.3) Vij = Q l VQ J , 

and Qj is the operator projecting into the j particles-j holes sector of the Q-space. 

The above corresponds to a partition of the Q-space into r sets, each one of these 
embodying a specific class of excited states: thus to the fc-th set are associated the fc 
particle-fc hole excitations. Actually, assuming the dominance of two-body forces in the 
residual interaction V (the operator connecting one set to the other), only excitations 
corresponding to an even number of particles (and holes) are actually filling the sets, as 
indicated in Fig. 1. The propagation of the system into the fc-th set is then described by 
the operator Gk, which, when k = r, reads 

(4.4) G r = ' 



whereas for k < r obeys the recurrence relation 

1 



(4.5) G k 



E - h kk - Vk,k+\Gk+iVk+\,k 
with 

(4-6) h kk = H QkQk + V QkP C^-^) V PQk . 
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Let us now, for convenience, pictorially describe a set as a box (see Fig. 1): then cq. (4.2) 
holds valid only if in the expansion the occupation of a box occurs when the previous 
box has also been occupied. 

Notice also that in a nucleus the expansion (4.1) is finite since the excited states 
with the greatest complexity correspond to the case r = A, A being the mass number. 

Since the mean field contains the average, the average of the corrections to the 
mean field vanishes. In addition because the processes occurring beyond the mean field 
are random, the average of the corrections to the shell model arising from any individual 
set vanishes as well. Therefore we should have 

(4.7) (AE) = (AEh) + (AE 2 ) + ■■■ + (AE A ) = 
and 

(4.8) (AE n ) = 0. 
However the average of the square of the corrections 

(4.9) ((AE) 2 ) = (J2^E n AE m ) 

has no reason to vanish. Furthermore, owing to the randomness of the Q-space physics, 
the nuclear matrix elements entering into the definition of AE n will also be random, 
hence the cancellation of the off-diagonal elements in (4.9), leaving 

(4.10) ((AE) 2 ) = <]T(A£„) 2 ). 

n 

To ensure the fulfillment of the requirement 

(4.11) (A£?i) = (AE 2 ) = ■■■ = (AE A ) = 
we redefine AEi according to 

(4.12) AEi = (^VoidVw^) - [(^iVkG^ol^)]^ , 

and likewise for the terms with r / 1. In (4.12) <jP A represents the nuclear mean field 
ground state wave function and the square brackets mean energy averaging. 

We are now in a position to calculate, according to the formula (3.21), the correc- 
tions to the mean field, provided the eigenfunctions ipk,a and the eigenvalues €k, a of the 
operator G^T 1 are known. We thus assume that the equation 

(4.13) G^Vfe.a = (E- e k , a )i> k , a 

can be solved for any value of the index k. 
Let us then first consider A£i: 



(4.14) ASl=E KV^W 



v l(y^|vWl 2 
v E - ^ 



AV 
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where e 2( g and \1p2p) are the eigenvalues and eigenvectors, respectively, corresponding to 
the 2p-2h excitations. 

Since in the present scheme only the square of the corrections to the mean field, 
namely 

(4-15) {SEhf =[{AE 1 )*] AV -[(AE 1 )] 2 AV , 

are meaningful, we are left with the task of evaluating 

<0 2/3 im>i 2 i(v>2 7 im)i 2 



(4.16) (SE.Y 



£ 

01 



(E - e 2[j )(E - e 27 ) 



AV 



£ 



Mia 

E — e 2 /3 



AV 



Introducing an average excitation energy e 2 for the 2p-2h states, the first term of the 
above equation can be recast as follows 



£ 

a,0 



(E - e 2a )(E - e 2/3 ) 



AV 



(4.17) = 



E-e 2 



£<<^L \V\M (ika I V\cj> A ){4>A \V\hp){^2p\V\<l>\) 
a,0 



AV 



The random phase average of the quantity in the square brackets is next obtained using 

(4.18) (AA*BB*) = (AA*)(BB*) + (A*B)(AB*) + quarticterms; 

quartic terms will be neglected. 

Now, the first term (AA*)(BB*) cancels exactly the second term in eq. (4.16). The 
evaluation of the second term in eq. (4.18) then yields 



(SE^ 



(4.19) 



E-e 2 
AEh 



E 



i\V\lp2 a )\ A v] 



(E - e 2 ) 2 D 2 



t [\{4> A \v\^)\ 



2 1 2 

AVI 



where AEi, essentially given by e, is the interval over which the energy average is carried 
out. The quantity D 2 is the energy distance between two neighboring 2p-2h states and, 
finally, -02 is meant as a representative of the 2p-2h excitations. 
We next consider AE 2 - 

AE 2 = (cj> A \V 01 G 1 V 1 2G2V2iG 1 V w \<p A ) 

y (^aj^M y>w\v\M ( ^ |y^o ) 

E — 69^ E — Car E — eo^' 



(E - l 2 ) 2 (E - e 4 ) 



E — 64,3 E — 62a' 

E^A^I^o)(02o|^|^> 

a,0 
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(4.20) 



■^ha\V\il>Afi){il>Afi\V\lh a ) 



J2\(4> A\v\^ a )\ 2 \(^ a \v\^)\' 

a, 



AV-i 



where we have set a = a', introduced an average energy £4 for the 4p-4h excitations 
and subtracted in the last line the average of AE 2 as written in the first line in order to 
achieve (AE 2 ) = 0. Focussing now on the square of the modified AE 2 = SE 2 we get 



(6E 2 f 



1 



1 



E &A\V\^ a )(i>2 a \V\<P° A ) 
a,p,a',f3' 



(e - i 2 y (e - e 4 ) 2 
^2 a \v\^){^m 2a ){ (4> A \v\^ 2a ,)^ 2a ,\v\cb A ) 



(MV\M(^\V\^ 2a ,)\ - J2 [\^A\VH2 a )\ 7 

' a,/3,a',/3' 



(4.21) 



|(^|V"|-02a')| 2 |(V'2a|V"|V'4/3)| 2 |('02a'|^l'04/3'}| 5 



AV 



AV 



To calculate the average we again invoke eq. (4.18). As before the first term cancels the 
subtracted term and one is left with 



(4.22) (6E 2 ) 



1 



1 



(e - i 2 y (e - e 4 ) 2 



Y,\^A\V\^ a )\ 4 \{^ a \V\^ f3 }\ 4 

a,0 



AV 



By comparing the above result with the one expressing (SE1) 2 , eq. (4.19), it follows 

1 1 



(4.23) 



(SE 2 y 



(E - e 2 )2 (E - e 4 ) 2 



(SEifK^V^)] 4 



from where an expression for the expansion parameter is emerging, the expansion being 
however a finite one as previously emphasized. Likewise ip 2 , in the above ip^ is describing 
a typical 4 particle-4 hole state. 

The generalization of eq. (4.23) to an expression connecting (5E n ) 2 and (5E n _i) 2 is 
then easily found to read 



(4.24) [{5E n f 



1 



1 



(E - e 2n ) 2 (E - e 2n _ 2 )- 



[(SE^f] [|^ 2 „- 2 |F|^„)| 4 ]. 



The results obtained up to this point have a general validity and can be applied to 
any system thus justifying the concept of mean field: indeed they have been deduced in 
the framework of the statistical theory. To further proceed a detailed description of the 
system one wishes to study is required: as an example we shall consider nuclear matter. 
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5. Explicit expressions for the projection operator and the HF theory 

We need now an explicit form for both the projection operator P and the NN 
interaction in order to obtain the true energy (E) and the mean field one (Eq) (the shell 
model energy for a nucleus). 

First we search for a suitable projection operator. Suppose, for this purpose, that 
for a given interaction the HF problem has been solved and let {0hf}> with the index 
i = 1,2,..., be the single particle wave functions of the HF orbitals, which form an 
orthonormal complete set for the Hilbcrt space of a one particle system. Out of the 
{<^hf) we can build an infinite set of Slater determinants {xhf! wn ich also form an 
orthonormal complete set for the Hilbert space of the nucleus. 

A natural choice for the projection operator is then 



M 

(5-1) ^ = £|Xhf)(XhfI 

i=l 

which, for M = 1, reduces to 

(5-2) ^HXhfXXhfI. 

|Xhf) being the ground state HF determinant of the nuclear system. In the following 
for sake both of simplicity and of illustration we shall stick to the case M = 1, although 
generalizations to larger values of M, while cumbersome, should be worth exploring. 

Notice that, with P given by (5.2), the wave functions (P^f) and turn out to 

be proportional to each other and to the ground state HF determinant |Xhf) as they 
must since the P space has only one member. 

With P given by (5.1) and the mean field by ii (eq. (3.12)) one gets 

(5.3) E = E HF + (xl i F\VQQV^\xl iF )- 



Eq — E — e 
where 

(5.4) E HF = (XhfI#pp|Xhf>- 

and 



E n -e-E 



(5.5) V = H Xl 

\ Eo-e-H QQ 

(see eq. (3.13)). Now writing writing Q = 1 — P in Eq. (5.3), one obtains 

(5.6) E = £ HF + ^^ {(XhfI^IXhf) - KXhfI^Xhf)! 2 }. 

Let us now specify the above general expression to the case of a Fermi gas (infinite 
nuclear matter), for which |xhf) = l-^ 1 )' 1^) being the wave function of a Fermi sphere 
of radius k F (the Fermi wavenumber). In infinite nuclear matter the quantity in curly 
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brackets in eq. (5.6), which we refer to as the variance of the residual interaction V, 
simplifies considerably. Indeed is (the states \n) labelling the spectrum of the Fermi 
sphere) 

(F\VV*\F) = Y^(F\V\n)(n\V*\F) 

n 

- (F\V\F)(F\V^\F) 
(5.7) + (F\V\l P -lh)(lp-lh\V^\F) 

lp-lh 

+ Y (F\V\2p-2h)(2p-2h\V t \F) 

2p-2h 

= \(F\V\F)\ 2 +(3 2 , 

the piece related to the lp-lh states essentially vanishing as required by the Brillouin 
theorem [7] and further terms in the right hand side being neglected because of the 
dominant two-body character of the residual interaction V. We indicate the variance of 
the residual interaction V, defined in eq. (5.6), with (3 2 to remind that in nuclear matter 
such a quantity turns out to be positive. To perform a reliable evaluation of [3 2 is of course 
quite hard because the residual effective interaction is an intricate operator: indeed, V is 
not only energy-dependent, but dependent upon the energy averaging procedure as well. 
In addition it is highly-non-linearly connected to the hamiltonian acting in the Q-space. 
However, even in abeyance of an explicit expression for the variance, the insertion 

(5.7) into (5.6) leads to the important result 

(5.8) Eq £* E RF 1 



Eq — E — e 

which embodies all the energies characterizing our problem, namely the HF, the mean 
field (which would correspond to the shell model in a finite nucleus) and the exact one. 

From (5.8) one sees that even if a residual interaction V is given and the associated 
HF energy £hf is calculated, still the energy E of the system cannot be obtained since 
the mean field energy Eq remains to be fixed. However, if as a first orientation one sets 
Eq ~ E, then the remarkable result 

(5.9) Eq = E HF - tL. 

follows, showing that the mean field energy is lower than the HF one, the parameter e 
being of course positive. This result holds only if \E -i?|<e. 



6. Occupancy of the mean field ground state 



The projected wave function P*f? is just a component of the full wave function 'J: 
how big a fraction of \P is contained in P 1 ^ is obviously a quantity of much interest. To 
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estimate it let (P^\P^>) = S 2 and search for an equation for S 2 in terms of the other 
quantities characterizing our problem. Since 

(6.1) S 2 = 1 - (Q*|Q*) 

using eq. (3.18) for (Q 1 ^) one obtains 



(6.2) 



S 2 = l-($o I V PQ ( fl,^ ) V QP \$o). 



Take now the Fermi gas model for i. e. \F). Then $o = S\F) and 
(6.3) 



S 2 = 1-S 2 (F\V PQ 



E-h 



QQ 



Vqp\F) 



which, with some algebra and neglecting the weak energy dependence of /iqq, can be 
cast into the form 



S 2 ^l + S 2 (F\V PQ — 



d 



dE \E — Hqq 



V QP \F) 



l + S 2 



(6.4) 

Because 
(6.5) 
and 
(6.6) 

it follows that 

S 2 = 1 
(6.7) +S 2 



-(F\ 



dE 
dVp Q 



V PQ 



1 



E - flQQ 



Vqp \F) 



dE E — IlQQ 

dV PQ _ 



Vqp\F) - (F\V PQ 



dV QP 



E-h QQ dE 



\F) 



dE 



Vpq 



1 



dV Q P 



dE^E-J^ V ^ 



2(E -e-E) 
1 

2(E -e-E) 



Eq — € — E 



(F\V PQ 



E - h QQ 



Vqp\F) 



or, by virtue of (3.21), that 

5 2 = 1 + 5 2 
(6.8) = 1 - S 2 



% {E - £o) + W^E 
dE 



+ ■ 



dE Eo-e-E, 

The explicit dependence of Eq upon E is provided by (5.8) which, when inverted, yields 



(6.9) 



-En — 



1 



E + e + £ HF ± V(E + e-E HF ) 2 +Ap 2 
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the minus sign in front of the square root being taken, since for vanishing residual inter- 
action (/3 2 = 0) the mean field should reduce to the HF field {En = Ehf)- Hence 



(6.10) 



dE 
dE 



E 



2dp 2 /dE 



11 -Ehf 
2 2 ^/(E HF -E-e) 2 + 4(3 2 ' 



which, inserted into (6.8), finally leads to 
5 2 = 



3 1 E HF -E-e- 2df3 2 /dE 
2 + 2 V /(S H 7"- E - e) 2 + 4/3 2= 



En — E - 



(6.11) 



3 1 E HF - E-e 
2 + 2 ^E^~E- e) 2 



2d(3 2 /dE + ^E -E HF 



-i -l 



4/3 2 



/? 2 



namely to the equation we were looking for. We thus see that the occupancy of the mean 
field ground state wave function S 2 is expressed in terms of the HF, mean field and exact 
energies, with an additional dependence upon the energy averaging parameter e, the 
residual interaction 1 and the derivative of the latter with respect to the energy. Since 
all these quantities are either explicitly evaluated (£hf, Eq and E) or fixed, in principle, 
both by the experiment (e) and by our theoretical framework (/3 2 ), eq. (6.11) provides an 
important check for the consistency of our approach, also because independent estimates 
on the depletion of the ground state wave function are available for lead, a nucleus for 
which nuclear matter is a reliable model, the present approach can be tested to some 
extent against the experiment. 

7. A simple model for nuclear matter 

To gain an insight into the effectiveness and self-consistency of the formalism de- 
scribed above we develop here a simple, schematic model. The model relies on two 
equations which, for sake of convenience, are here repeated. The first of these is eq. (5.8) 



Eq = -Ehf + 



3 2 



En — E 



The second is obtained from Eq. (4.19), i. e. 

2 



(SE 1 



1 



E-e 2 



which we heuristically approximate (admittedly roughly) by 



(7.1) 



/3 4 



(E - aef 



where a is a parameter which estimates the average energy of the 2p-2h states involved 
(we set 62 = «e for convenience) and importantly helps to correct for the poor approxi- 
mation of /3 4 for the sum over the 2p-2h states 7 in eq. (4.19). Taking the square root, 
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from eq. (7.1) it follows 



(7.2) 



E — Eq — dz 



E - at 



If now one combines eq. (5.8) and (7.2), using the upper sign (plus) one obtains 
a lower bound Ei for the energy of nuclear matter, while using the lower sign (minus) 
one obtains an upper bound E u . Clearly, one gets also differing values of the mean field 
energy to be denoted by Eq and Eq, respectively. We ask then whether there is a choice 
of e, a and [3 such that 

i) the two mean field energies E and Eq turn out to be the same over a range of 
densities (or of Fermi momenta kp) of significance for nuclear matter; 

ii) the "experimental" values of the binding energy, saturation density and compres- 
sion modulus of nuclear matter, namely B.E./A = —16 MeV, Uf — 1-36 fm" 1 and 
K = 23 or 16 MeV (the two values refer to a hard and a soft equation of state, 
respectively) [8] are accounted for in a sense to be later specified; 

iii) the spectroscopic factor given by eq. (6.11) turns out to be less than one both 
on the lower bound (where its value Si is associated with Ei) and on the upper 
one (where its value S u is associated with E u ) and furthermore not too much at 
variance with existing estimates. 

To accomplish this program we proceed by first eliminating Eq and E Q from eq. (5.8) 
and (7.2). This yields an equation for (3f and /? 2 : 



(7-4) 01 = -^j^ { [£hf - e(a - 1)] - VI^hf - e(a - l)] 2 + Ae{E u - £ HF )} , 



for the upper bound. Note that the solutions with the plus sign in front of the square 
root have been discarded, since they lead to an incorrect limit as e — ► 0. The energy 
derivatives of (7.3) and (7.4)are then easily obtained, namely 





dE 



(7.5) 



for the lower bound, and 




2 {[#hf - e(a - 1)] - ^[E KF - e{a - l)] 2 + Ae(E u - E HF )} 



(7.6) -(E-ae) 



V^hf - e(a - l)] 2 + Ae(E u - £ HF ) 
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for the upper bound, These, when inserted into (6.11), provide an expression for the 
spectroscopic factor, associated with the lower and the upper bound, respectively, in 
terms of Ei (E u ), E l (Eft), Eu F , a and e. 

From eqs. (7.3) and (7.4) we note that should appropriate values for Ei (or E u ) and 
-Ehf be available one would then obtain a value of fif and [3 2 for given values of a and 
e. These values of [3 2 can in turn be used in eq. (5.8) to obtain E and Eft. We agree 
that a solution of eq. (5.8) and (7.2) occurs when these upper and lower bounds for the 
mean field are equal. 

The problem remains of how to choose meaningful values for Ei and E u and how to 
fix Ehp. The first of these questions is answered by taking 



(7.7a) Ei=S- W/2 

(7.7b) E U =S + W/2, 

where the energy 

(7.8) £= [-16 + 39. 5(fc F - 1.36) 2 ]MeV 



incorporates the present knowledge on the ground state energy of nuclear matter, quoted 
in ii), as extrapolated from finite nuclei and W should be viewed as the fluctuation energy. 

The value of E^p as a function of Uf is obtained from an assumed, schematic two- 
body interaction (Appendix A). The resulting Hartrcc-Fock energy for a Fermi gas is 
shown in Fig. 2 and, for a reduced scale in kp, in Fig. 4. From Fig. 2 we see that for 
our simple two-body force a minimum of the binding energy as a function of &f occurs 
at Uf = 1-78 fin -1 , the corresponding energy per particle being —7.3 MeV. 

The potential can also be used to calculate the bare value of fi 2 (Appendix B). 
Obtained values are shown in Tables I, II and III: they range from 4.5 x 10 4 to 7.0 x 10 4 
MeV 2 /nucleon. 

In discussing our findings it help first to observe that actually, for a given fluctuation 
energy W, a whole set of values for the parameters a and e exist such to satisfy the 
requirement i), namely E — Eft, at Uf = 1-36 fm _1 . They lie on the curves displayed 
in Fig. 3. Furthermore a and e should also be such to fulfill the constraint 

(7.9) e 2 = at > 20MeV, 

which represents a fair estimate of the lower limit for the excitation energy of the 2p-2h 
states. As a consequence the acceptable values for a and e are then restricted to the 
domain to the right of the dotted line in Fig. 3. 

In this region a and e should be selected in such a way to comply with the require- 
ments ii) and iii). This turns out to be possible and our results are shown in Fig. 4 and 
Tables I, II and III. 

However there one sees that as one moves away from kp = 1.36 fm _1 the equality 
between Eq and Eft is no longer exactly satisfied, although the two mean fields remain 
rather close to each other at least in the range 0.9 < kp < 1.5 fm _1 providing the values 
of W are "moderate" . Indeed, quantitatively, requirement i) can be reasonably satisfied 
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in the above quoted range of Uf for W < 4 MeV only. For larger fluctuation energies E l 
and Eq differ too much: an orientation on the size of the error around the mean field in 
thus obtained. 

Concerning the energy averaging parameter, e, it turns out to be about f MeV larger 
that W for all the cases listed in the tables. The values of the spectroscopic factors Si and 
S u are quite stable, while the values of the effective interaction f3f and (3^ grow rapidly 
with increasing W . Most importantly, their values differ by three orders of magnitude 
from the bare (3 2 . This occurrence partly stems from the renormalization (see eqs. (3.f3) 
and (3.14)) induced by the energy averaging. The random phase averaging also gives rise 
to a reduction in the magnitude of the residual interaction. Calculations, while rough, 
still indicate that these two effects appear indeed sufficient to produce the observed sharp 
reduction. 

8. Closing comments 

The statistical theory of the mean field presented in sections 1-4 is based on two 
propositions, ft is assumed that the mean field is the slowly varying component of the 
nuclear interaction, which can be obtained by taking an appropriate energy average. 
Secondly it is suggested that the matrix elements of the residual interaction are random, 
so that their average is zero. A formalism incorporating these ideas, borrowed from 
statistical reaction theory, was developed and explicit expressions for the mean field and 
the "fluctuation" away from the mean field was obtained. 

Next an expansion of the fluctuation energy in terms of increasing excitation com- 
plexity leads, after averaging, to formulas for the corresponding contributions to fluctu- 
ation energy. Worth emphasizing is the novelty of this approach: indeed the corrections 
to the mean field ground state energy of the system rather than being ordered according 
to the power of a coupling constant or according to the number of hole lines occurring 
in a given diagram, as it is done in standard theories, are here organized in terms of the 
complexity of the Q-space states reached via the residual interaction acting among the 
constituents of the system. 

The presentation of this approach was followed in sect. 5, 6 and 7 by a simplified 
version and by a schematic model. The most remarkable result emerging from the anal- 
ysis of the model has been the sharp reduction of the effective strength of the residual 
interaction. But in addition is the overall reasonableness of our results which is quite 
encouraging. 

Of course much remains to be done. The quantitative connection with the underlying 
nuclear forces has not been exhibited. The evaluation of the matrix elements for finite 
nuclei was not carried out and a better understanding of the energy average needs to be 
achieved. We have so far only a schematic model. What is needed is a complete and 
thorough evaluation which, on the basis of the obtained results can be expected with 
confidence to be successful. 

Application to excited states is also indicated. In this case the smoothing function 
used in reaction theory can be used instead of the one of section 4. This would lead to 
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a complex mean field and the excited states would have a width corresponding to the 
probability of the splitting of the state by the residual interaction. The width would then 
measure the extent of the splitting. 

Of crucial importance is also to test the rate of convergence of the complexity ex- 
pansion. Three elements are expected to favour in general a fast convergence rate for the 
latter: 

i) the energy averaging parameter e, if it is large enough; 

ii) the subtraction of the average interaction from the bare one; 

iii) the reduced overlap between the ground state and the complex wave functions 
belonging to the Q-space. 

A detailed investigation will be invaluable in sheding light on this issue. In this 
connection calculations are presently in progress. 



Appendix A. 



The bare interaction and the HF theory 

To implement the program outlined in Section 7 we need a NN interaction to fix the 
HF energy of nuclear matter. For illustrative purposes we choose the following simple 
NN interaction 



(A.l) 



V(r)=g A - 



9b- 



e -W I + p x 



which embodies a short-range repulsion in the first term and an intermediate range 
attraction in the second one (all the parameters are positive and we require \ia > Ms)- 
The latter is taken of Majorana type, hence the occurrence of the exchange operator 



(A.2) 



P, = 



1 + Ti • T2 1 + <?i ■ <?2 



built out of the spin a and isospin f operators. Thus (A.l) contains the main features 
needed to account for the saturation of the nuclear forces, reflected in the existence of a 
minimum in the binding energy per particle (B.E. /A) versus kp curve, but for the tensor 
force that here, for the sake of simplicity, is neglected. 

The HF energy for the interaction (A.l) is easily worked out and leads to the fol- 
lowing expression for the binding energy per particle [9] 



(A.3) 



B.E. 



3 h 2 k 2 F 
5 2m 
1 



8Trk f 



1 

+ - 

7T 



9aA 



9 a 37T g B \ 3k F 



47T 



:9b ^ B 



gAl^A arctan 



2fc F \ 3 

+ t;9b^b arctan 

Ma / 2, 



2 9B 



2k F 
Ms 
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8nkF 



9AI& 



3 + 



4fc F 



log 1 



( 3 + it) log ( 1 + 7f 



2k%/3n 2 . 



where g — ^ F / 

The parameters characterizing the potential (A.l) might be fixed, for example, by 
accounting for the "experimental" nuclear matter values previously quoted. One succeeds 
in doing so with the following choice 



(A.4o) 
(A.46) 



9A 



3.43fm , 
2460McV fm, 



fiB = 1.63fm 
g B = 898MeV fm. 



Worth noticing is that the range of the repulsion obtained with the fitting procedure 
is rather close to the one associated with the exchange of a ui meson (3.97 fm -1 ), whereas 
the range of the attraction turns out to be intermediate to the one arising from the 
exchange of a pion and of a a meson (0.71 fm -1 and 2.79 fm -1 , respectively). 

Here, we rather prefer to choose the parameters in such a way to have too little bind- 
ing energy at too large a density in the HF frame, in order to conform to a shortcoming 
common to many nuclear matter calculations, the purpose being to ascertain whether 
the present theory is capable to improve upon the HF results. Of course, there exists a 
variety of ways for reaching this scope: In view of the rather realistic values of the ranges 
Ha and /j,b (see (A. 4a)) we change the coupling constants. We thus take, as a rather 
extreme example, 

(A.5) g A = 740McV fm, g B = 337McV fm, 

which, together with the values for \xa and \i B given in (A. 4a), yields a minimum of 
— 7.30 MeV for the binding energy at k F = 1-78 fm -1 , as it can be seen in Fig. 2, where 
the HF energy is displayed as a function of the density. 



Appendix B. 



Vacuum — > 2p 2h matrix element 

In the language of second quantization the matrix element we have to calculate reads 

(B.6)/3 2 =]T ]T (kMV\k+q,k-${k+qM-m\ki,k 2 ), 

spin fci,/c2<^F 
isospin |fe 1 +£|,|fc 2 _g1>fc j: . 

which, with standard manipulations, can be transformed into 

xB(k F -fc 2 ) \ 32g 2 2 1 +12g% 2 - -24g A g B - 



K + ^(A + 1 2 ) 2 A + Q^b + Q 2 
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(B.7) -8g 2 A 



+ 12g_ 



B~~2 



Ml + <? 2 /i| + |fci-fc 2 +QT 



-l2g A g B 



+ 



n + q 2 mI + |fci -k 2 + q\ 2 mI + q 2 p? A + |fci - fc 2 + gp 



(the matrix elements in eq. (B.6) are understood to be antisymmetrized) . The above 9- 
dimcnsional integral, by appropriate transformations, can be reduced to a combination of 
2- and 4-dimensional integrals, which can be numerically evaluated, yielding the results 
quoted in Tables I— III. These values turn out to be helpful in performing a comparison 
with later findings. 

The direct contributions, namely the first three terms in the right hand side of (B.7), 
are analytically evaluated. Considering, e. g., the third one, one gets 



A ir q J WFWF m O* 1 + ^ ~ kp ) 9 ~ ^ ~ kp ) 9 {kp ~ kl) 9 {kp ~ k2) 



(B.8) 



x 



/4 + q 2 mI + q 2 



A 

9 17 



A-r, — \ -=o (ft A arctan /i A - p, B arctan /2 B ) 

4fc F 1 Ma - MS 



+ 



i 



4 #4 - Al 



Ms(3 + Ms) arctan 



1 

Ms 



Ma(3 + Ma) arctan 



1 

Ma 



where ma,b = Ma.b/2/cf. 

The other terms are obtained with obvious substitutions. 

The exchange contribution is somewhat more involved. As an example we consider 
(again the other pieces are obtained through obvious interchanges) 



irg 



/dk\ dk 2 
(2tt)3 (2tt) 



■dg 



© (|ifei + gl - **■) e (|fc 2 - q\ - k F ) e (** - h) e (fc F - k 2 ) 



(B.9) 



(2 



iMA + 9 2 M| + |fci -fe + g? J 

4 9 00 /-oo 2 /■ 



drre 



fJ-Br ■ 



Z=0 

k F 



dkk z ji(kr) I dxPi(x)Q(q 



Ma + 9 

2 1 k 2 + 2qkx - k F ) 



jo(qr) 



where a partial wave expansion, in terms of the Legendre polynomials Pi and the spherical 
Bessel functions ji, has been performed. The above integral, for q > 2kp is easily 
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calculated yielding 



2 f°° e~^ BT ( y 

P exc = A^- / dr - — (sm(kFr) — kpr cos(kFr)) 

n Q Jo r v ' 



(B.10) x f 

J% 



Jo 

X / d Q 2 q , 2 sin (g r )' 



The remaining integrals are numerically evaluated. 
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Table I. - In the second row the value of Eo as a function of k F is reported; in the third 
row one finds the bare value of f3 2 , as given by eg. (B.7), while in the fourth and fifth rows 
the estimated renormalized values on the lower and upper bounds, respectively; in the last two 
rows, the corresponding values of the spectroscopic factor are displayed. The values in this table 
correspond to W = 2 MeV, e = 3 MeV and e 2 = 22 MeV. 



k F (fni -1 ) 


1.2 


1.36 




1.5 


E (MeV) 


-15.3-^-15.5 


-16.3 




-15.6^-15.2 


p 2 

Pi (MeV /nucleon) 

Pi 


4.5-10 4 
27.0 
54.2 


5.2-10 4 
26.3 
49.0 




7.0-10 4 
22.8 
36.4 


Si 

Su 


0.32 
0.64 


0.32 
0.62 




0.30 
0.58 


Table II. — As in table I: 
€2 = 20 MeV. 


The values in this table correspond to W = 


3 MeV, e 


= 4 MeV and 


k F (fin -1 ) 


1.2 


1.36 




1.5 


E (MeV) 


-15.5-^-15.9 


-16.6 




-15.9-^-15.2 


Pf (MeV 2 /nucleon) 

Pi 


4.5-10 4 
36.1 
79.6 


5.2-10 4 
35.3 
69.9 




7.0-10 4 
30.7 
49.8 


Si 

S u 


0.34 
0.69 


0.33 
0.67 




0.32 
0.63 


Table III. - As in table I: 
e 2 = 19 MeV. 


The values in this table correspond to W = 


4 MeV, e 


= 5 MeV and 


k F (far 1 ) 


1.2 


1.36 




1.5 


E (MeV) 


-15.7^-16.4 


-16.8 




-16.1-^-15.2 


P 2 

Pf (MeV 2 /nucleon) 

Pi 


4.5-10 4 
45.7 
108.6 


5.2-10 4 
44.7 
91.9 




7.0-10 4 
39.2 
62.7 


Si 
S u 


0.35 
0.73 


0.34 
0.72 




0.33 
0.66 
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Qo 


> 


Q 2 


-> 


Q 4 


> 


< 


< 


< 









g. s. 2p-2h 4p-4h 

Fig. 1. - The partition of the Hilbert space of nuclear matter in sets of increasing complexity. 
The first box on the left defines the P-space, the second one embodies the simplest states in the 
Q-space and so on. 




0.5 1 1.5 2 2.5 

k F (frn -1 ) 



Fig. 2. - The binding energy per particle in the HF approximation (formula (A. 3) of the text) 
for the potential (A.l) and the parameters given by (A. 5). 
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Fig. 3. - The loci corresponding to E l = Eq in the plane (a, e), for W — 2, 3 and 4 MeV. Also 
shown (dotted line) is the curve along which the average energy of the 2p-2h excitations is 20 
MeV 
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1 1.2 1.4 1.6 



k F (frrT 1 ) 



Fig. 4. - The binding energy per particle in the HF approximation (dot) and in the present 
approach: The solid lines represent the lower and upper bounds of fluctuations for the energy 
E, whereas the dashed lines give the corresponding mean field energies, (a): W = 2 MeV, e = 3 



